Anisotropic current-induced spin accumulation 
in the two-dimensional electron gas with spin-orbit coupling 
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We investigate the magnetoelectric (or inverse spin-galvanic) effect in the two dimensional electron 
gases with both Rashba and Dresselhaus spin-orbit coupling using an exact solution of the Boltzmann 
equation for electron spin and momentum. The spin response to an in-plane electric field turns out 
to be highly anisotropic while the usual charge conductivity remains isotropic, contrary to earlier 
statements. 



I. INTRODUCTION 

Spin phenomena in semiconductor structures lie at the 
very heart of the emerging field of spintronics and are 
to a major and still growing direction of solid-state re- 
search. Among the plethora of concepts and ideas, the 
magneto-electric effect in the two-dimensional electron 
gas (2DEG), sometimes also referred to as the inverse 
spin-galvanic effect, has attracted particular interest 
from an both experimentali'^'^'^'^ and theoretical^i^iSiiSiiiS 
point of view. This effect amounts in spin accumulation 
as a response to an applied in-plane electric field and is 
therefore a possible key ingredient towards all-electrical 
spin control in semiconductor structures, a major goals 
in today's spintronics research. 

We find it important to emphasize that the phe- 
nomenon studied below is not the intrinsic spin Hall 
effectiiii^ though it would manifest in a very similar man- 
ner. The mechanisms involved in these two phenomena 
are, however, very different. In particular, the momen- 
tum relaxation due to the impurity scattering is necessary 
for the spin accumulation investigated in present paper, 
contrary to the intrinsic spin Hall effect where ballistic 
regime is preferre d^^i^^ . 

In this paper we investigate the magnetoelectric effect 
in 2DEG's with both Rashba and Dresselhaus spin-orbit 
coupling. Our study is based on an exact analytical solu- 
tion to the Boltzmann equation for electron spin and mo- 
mentum in the presence of (5-function shape static impu- 
rities. Regarding the spin degree of freedom, our solution 
neglects off-diagonal elements of the semiclassical distri- 
bution matrix in the eigenbasis (or helicity basis) of the 
single-particle Hamiltonian, an approximation which is 
shown to be valid at sufficiently high temperatures com- 
mon in experiments. 

As a result, the magnetoelectrical effect, i.e. the mag- 
netic response to an in-plane electric field, turns out to be 
highly anisotropic, whereas the usual charge anisotropy 
remains isotropic, contrary to earlier statementsiS,. 

The paper is organized as follows. First, we present 
an analytical solution of the semiclassical spin-incoherent 
Boltzmann equation. Second, we solve the spin-coherent 
kinetic equation. Finally, we discuss the plausibility of 
solutions obtained and apply them for investigations of 
the electric current and spin accumulation. 



II. SOLUTION OF THE KINETIC EQUATION 

We consider the Hamiltonian as a sum of the kinetic 
energy and two spin-orbit coupling terms: RashbaM and 
Dresselhaus^^. Then, the Hamiltonian takes the form 



H = 



2m 



+a{axky-aykx)+f3{(7xkx-(7yky). (1) 



Here, ax^y are the Pauli matrices, k^^y are the electron 
wave vectors, and m is the effective electron mass. Intro- 
ducing the angle 7^ so, that 

akx 



tan7fe 



f3ky 
Pkx + aky ' 



(2) 



we obtain the following spinors as the eigen functions of 
the Hamiltonian H 



ik^x+ikyy 
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The energy spectrum has the form 



E. 



± 



2m 



± \J {akx + PkyY + {[ikx + akyY 



(3) 



(4) 



where k = \ kl 



fcy . The velocity matrix in the helicity 
basis ^ is not diagonal, and its elements read 

-y^"^'^^-* — hkx/m ± (/3 cos + Oi sin jk)/fi, (5) 
w^2(2i) ^ ±2(/3sin7fc - acos7fc)/n; (6) 
11(22) ^ fiky/m± (acos7fe + I3s\\\^k)/K (7) 



12(21) = ±i(Q,sin7fc - /3cos7fc)/?i. 



(8) 



For the diagonal elements of the velocity matrix we also 
use simplified notations v-'--'-'-^'^^ = v±. 

In the following, it will be convenient to use polar co- 
ordinates kx — kcosO, ky — ksinO. Then, the spectrum 
read E± = ?i^fc^/(2m)± \ k \ ng, where 



Kg = y/a^ 



■ -f 4a/3 sin 6 cos 9 



is the generalized spin-orbit interaction constant for a 
given direction of motion. The wave vectors for a given 
energy read 



2mE 



(9) 
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and the expression for 7^ takes the form 
OL cos + /3 sin G 



tan7fc 



(i cos 61 + a sin d 



We model the influence of impurities using (5-potential 
scattering V = Cd{x,y). The scattering probability be- 
tween the states with fc, s and k' , s' then read 

w{ks; k's') = [1 + ss' cos(7fe - ^k')] 5{Esk - Es>k'), 

rriT 

(10) 

where r is the relaxation time which relates to ^ as 
T = h^/{mO. Note, that the spin-dependent factor 
in Eq. pH)) does not depend on the particular form of 
V{x,y) in the Born approximation. We expect, there- 
fore, that such a simple scattering model as (5-function 
shape static impurities is quite reliable for the descrip- 
tion of spin-dependent transport phenomena. 



A. Spin-incoherent kinetic equation 



In general, the equation has the form 
9/s(k) 



dt 



V, ■'r + (-eE) ^ ' 



dr 



hdk 



dfs 
dt 



(11) 



Here, s = ± is the spin index, and Vg are the diagonal 
elements ([5]) and (O of the velocity matrices. 

To solve Eq. pT|) we follow the standard procedure 
widely spread in literature (see Ref.-^, or any textbook 
on solid state physics). We write down the distribution 
function as fs = fs + fl + : where /° is the Fermi 
distribution and /^'^ ^ Then, at zero temperature 
gradient and constant electric field in the linear responce 
regime Eq. (fTT|) takes the form 



eEv.s 



dEsk 



= St[/s(k)]. 



(12) 



Assuming elastic scattering fullfilling the microreversibil- 
ity condition the scattering operator can be written as 



St[/.(k)]=^ j 



<Pk' 
(2^)2 

/l'(k' 



{«;(ks;k's')[/l(k)+/2(k) 
l-/?(k')]}, 



(13) 



where ^(ks; k's') is given by Eq. (fTT 

We would like to emphasize that besides abovemen- 
tioned conventional assumptions our solution is exact. It 
means that in contrast to a number of previous works we 
do not make any further approximation from this point. 

The solution of Eq. (I12|) has the form 



fl 

fs = 



reEvs 
TeE 



dp{Esk) 



dEsk 
df{Esk) 



(14) 



dE. 



sk 



(15) 



X [{ax cos 9 + bx sin 9)ex + {ay cos 6 + by sin 6)ey] 



Here, f} is the conventional solution of the Boltz- 
mann equation, whereas an addition is usually dis- 
carded since it is much smaller than f}. However, this 
anisotropic addition plays an important role in the spin 
accumulation mechanism as it is shown below. 

In order to define the unknown coefficients ax,y, bx,y, 
we substitute f}'^ into Eq. (fT^ (see Appendix for details) 
and get the following equations 



a 



13' 



a(3 



1 + 



4a/3 
{a' +13'? 



_ ^ _ 
2 



a 



(3' 



2a/3 



0' 



(16) 
(17) 



The equations for ay and by can be obtained from 
Eqs. p6|) - p?)) by the substitution ax — + by, bx — > ay. 
From these equations one can easily find that ax — by — 
-(a^ -I- /3^) and ay ^ bx = -2a/3. 

Finally, we substitute dS]), (O into Eq. ([14]), sum up f} 
and /2, and write down the solution of Eq. (fT2|) in the 
following elegant form 



/,^/° + (-eE)k^ 
m 



df{Esk) 
dEsk 



(18) 



This solution should be compared with the expression in 
the absence of spin-orbit coupling which can be formu- 
lated in terms of the velocity as follows 



/, = /o + (-eE)v,r 



df{Esk) 



dE. 



sk 



(19) 



where = Kk/m is the velocity. Note that this latter 
relation holds only in the absence of spin-orbit coupling. 



B. Spin-coherent kinetic equation 

Since we deal with the constant electric field only, the 
equation readai^. 



(-eE) 



dm 

hdk 



HJ{k) 




(20) 



Following Ref.^^, we rewrite Eq. ((20)) in the helicity ba- 
sis where the Hamiltonian is diagonal, and the equation 
takes the form 



f hi 


fl2 \ 


ieE djk 


f f21 


U21 


f22 J ' 


2h dk 


[hi 







fl2{E+ - 


E-) 


f2l{E 


--E+ 








where 



= St[/(k)], (21) 



d-fk . 

—— = — —2 — (e^; sm - By cos 0) . (22) 
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The collision term readsii 



X [k:(^, (fc) + K^^ f,,,, (k)]}, (23) 



where 



" irriT 



(24) 



Eq. p4)) was derived for 5-potential scattering introduced 
above, and each K^f^, contains either sin(7fe — 7^/) or 

1 ± cos(7fc — 7fc') depending on the product of ss'sis[. 



In order to simplify Eq. pi B[) . we assume kng ^ T, 
where T is the temperature. This assumption will be 
discussed later in details. Now, we expand the Fermi 
distribution function f'^{Esk) in terms of kng/T so, that 
Eq. (jll B[) in the linear responce regime takes the form 



eE 21 
2 ^ 



2 ^ 



\f21iE- - 



= Si[/(k)]. 

m 



Though Eq. ([^5]) is still quite cumbersome, it is easy 
to check that /12 = /21 = and /ii(22) given by 
Eq. represent the solution. The details of calcu- 
lations can be found in Appendix. Thus, the solution for 
the spin-incoherent Boltzmann equation is the same as 
for the spin-coherent kinetic equation at high tempera- 
tures. This can be explained in what follows. 



III. RESULTS AND DISCUSSION 

Note, that the off-diagonal elements of the distribu- 
tion function have essentially quantum mechanical ori- 
gin since they correspond to the off-diagonal elements of 
the density matrix. Classically, an electron can be in 
only one state of two, and, therefore, off-diagonal terms 
vanish here. In contrast, in the spin-coherent kinetic 
equation the off-diagonal terms can be essential. How- 
ever, in the real samples the quantum effects are neg- 
ligible at room temperatures because of the tempera- 
ture smearing. Indeed, the spin-orbit coupling constant 
Ke is of the order of 10~^^ eV • m for typical InAs sam- 
ples. To be specific, let us take n-type InAs quantum 
well containing the 2DEG used for photocurrent mea- 
surements at room temperature^. The parameters are as 
follows: a/ (3 = 2.15, mobihty is about 2 • 10^ cm^/ (V • s) 
and free carrier density is 1.3 • 10^^ cm^^. The latter 
allows us to estimate characteristic Fermi wave vector 
kp = V^Trrie ~ 3 • lO^cm"^. Thus, the spin-orbit split- 
ting energy kpKg is about 3meV, that is much smaller 
than Tioom = 25 meV, and our solution is suitable for 
description of a large variety of experiments. 

To study the spin accumulation we calculate the net 



spin density, whose x, y-components read 

where Sx{k,s) = |cos7fc, Sy = — |sin7fc are the spin 
expectation values. (We do not consider S'^-component 
since it is zero.) 

The integral over k can be taken easily making the sub- 
stitution e — E{s,k) and assuming that —df'^{e)/de ~ 
5{Ef — s). This assumption is reasonable with respect 
to the system studied in Ref.^* from which one can de- 
duce a Fermi energy of the order of 100 meV, which is 
clearly larger than room temperature. Since our solution 
is valid for temperatures much higher than the spin-orbit 
splitting energy kpKg, the inequality describing the ap- 
plicability of our results obtained below reads 

kpKe <^T ^Ep- 

The rest integrals over the polar angle can be taken 
analytically. After some algebra we have 

Then, the magnitude of the spin accumulation (S) — 
\/{Sxy + (Sy)^ is given by 

p T'jTnT I - — - 
(5) = ^^V"'+/'' + 2«/3sin(2Ee,). (28) 

It is interesting to note, that {S) depends on the direction 
of the electric field (see Fig. [ij, i. e. the spin accumu- 
lation is anisotropic. To our knowledge, this interest- 
ing feature was not noticed in the literature so far. If 
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FIG. 1: Spin accumulation (in arbitrary units) vs. direction 
of the electric field in polar coordinates for different Rashba 
and Dresselhaus constants: A — a — 3/3, B — q = 2.15/3, C 
— a = (3. The curve B corresponds to the real situation in 
n-type InAs^^. 



13 = {) then (5) = eEmTa/{2TTh^) that is in agreement 
with Refs.^'^. If the electric field is applied along the x- 
direction (the case studied in Ref.^), then the spin density 
is (S) — eEniT^a^ + fS^ I[2'k^). This result contradicts 
to Ref.^, where the spin density has some strange kinks 
as a function of mcP' jY? and m0^ /h^. However, one can 
see from Eq. (|28|) that the kinks in the dependencies of 



{S) on a (or /3) could not take place. 

Relying on the effect depicted in Fig. [1] the follow- 
ing novel spintronic device can be proposed. Let us 
attach two pairs of contacts to the 2DEG so that the 
first pair provides the electric current along the crystallo- 
graphic axis corresponding to the minimal spin accumu- 
lation (i. e. [110]-axis for [001]-grown InAs samplesl?), 
and the second one is connected along the perpendicular 
axis. Then, the spin accumulation depends strongly on 
which contacts the transport voltage is applied, and its 
anisotropic contribution can be extracted easily using op- 
tical methods^'"=ii^ or just measuring the magnetization. 
To give an example, applying the electric field of 20 V/ni 
(which corresponds to the current density of ImA/cm) 
we obtain the magnetization difference of the order of 
lO^/is per cm^, where is the Bohr magneton. This is 
comparable to the Pauli magnetization at the magnetic 
fields of a few gauss. It is also interesting to note the 
small characteristic switching time t ~ lO^^'^s of the 
device proposed. Therefore, besides the fundamental im- 
portance of such an experiment, our four terminal device 
could find some applications as a high-speed spin switch. 

At the end of the discussion, let us turn to the charge 
current whose density can be found as 



Vs/s(k). 



(29) 



ner as in p6)) . and after some algebra the conductivity 
tensor a takes the form 



'yy(xx) 



eV 



If 



(30) 



and, most surprisingly, <Jxy(yx) ~ 0. It is convenient to 
express <Jyy(xx) via the electron concentration pip in a 
2DEG with spin-orbit interactions 



mEp ( m 



f3' 



(31) 



Then, the conductivity takes much simpler form, namely 
^yy(xx) = e^nf-T/ra. This is the Drude formula, i. e. 
the conductivity is just a number (not a tensor), though 
the distribution function is anisotropic. This result 
is in contrast with the findings of Refii^, where a sub- 
tle approximation generated off-diagonal elements in the 
conductivity tensor—. This artifact is absent in our truly 
exact solution: The electrical conductivity is isotropic for 
any relation between a and /3. The latter simplifies es- 
sentially the relation between the spin and charge current 
densities, which takes the form 



P a 



(32) 



We find it useful to write down Eq. (|32p in the basis 
-^(1;1), -^(1;— 1). Then the relation between (S) and 
j takes the simpler form 



27r?i3e 





/3-f a 



(3- a 




(33) 



The integral over k in (j29p can be taken in the same man- 



From this equation one can see easily, that the spin accu- 
mulation is strongly anisotropic if the constants a and (3 
are close to each other. The reason of such an anisotropy 
is the angular dependence of the dispersion law In- 
deed, the spin-orbit splitting is different for different di- 
rection of the momentum. Therefore, the spin precession 
frequency depends essentially on the direction of the elec- 
tron motion. This leads to the anisotropic spin relaxation 
times (see e. g. Refi^), and, thus, the anisotropy of the 
spin accumulation occurs. In particular, the electric cur- 
rent of arbitrary strength applied along the [110] crystal- 
lographic axis does not lead to any spin accumulation at 
a = 13. The latter is due to the vanishing spin splitting 
along the [110]-axis (see e. g. Ref^^). 

Eqs. ((32|l and ((33|l are the main result of our work 
which can be applied directly to experimental studies of 
current-induced spin accumulation in [001]-grown InAs 
samples. The theory developed here is, of course, appli- 
cable to arbitrary oriented heterostructures after a few 
minor changes regarding the spin-orbit coupling terms in 
the Hamiltonian. 



IV. CONCLUSIONS 

In summary, we solved the semiclassical Boltzmann 
equation analytically for 2DEGs with arbitrary large 



5 



spin-orbit interactions of both Rashba and Dresselhaus 
type. Moreover, we demonstrated that this solution is 
also suitable for the spin-coherent case at sufficiently 
high temperatures common to experiments. Using this 
solution, we discovered the anisotropy of the current- 
induced spin accumulation, though the conductivity re- 



mains isotropic. Finally, our analytical study is expected 
to be a reliable starting point for further investigations 
of spin dependent electron transport. 

We thank S.D. Ganichev and D. Loss for useful discus- 
sions and acknowledge financial support from Collabora- 
tive Research Center 689. 



APPENDIX A: THE COLLISION INTEGRAL CALCULATIONS 

To take the integrals of the form / -^^^ F {k' , E gi k') it is convenient to perform the substitution Es'k' = where 
F(/c', Es'k') is a given function, and Es>k' is the dispersion relation 



1. Spin-incoherent case 

The first term of St[/s(k)] given by Eq. ((nX)) reads 

E / |^-(k^;k'^')[/l(k) + /.^(k)] = [/l(k) + /,^(k)]X:/ ^ J de'^ 



5{Es 



1 + ss 



, [a^ + 13^) cos{0 - e') + 2al3sin{9 + 9') 



KeKo' 



-[/i(k) + /2(k)]. 



2me' 



(Al) 



The rest terms of St[/s(k)] containing fl and given 
by Eqs. (fHl) and respectively read 



E 



dpk' 



E—seEx.y 
2hK0 



df(Esk) 



cos f 
sin^ 



dEsk 
2af3 



smt 
cos( 



(a2 + /?2)2 /|a2_^2 I 



4a2/32 



(A2) 



E 



df{E,k) 



EseEx^y 
2hKe 

x,y 



COS f 



dEsk 



2a(3 



sm6 bx.y + a 



x,y 



a" + /g^- I - I 
2^ 



(A3) 



Substituting Eqs. (jAl[) - (jA3p into the master equation 
(fT2)) one can easily establish Eqs. p6)) - (|17p for unknown 
coefficients ax,y, bx,y 



2. Spin-coherent case 

To prove the solution of Eq. (|25p we take the integrals 
in its right hand side in a similar way as before. Indeed, 
after the substitution of /ii(22) given by Eq. (fT8|) and 
/i2 = /21 = into the collision term we have 



St[/(k)] 11(22) 



/: 



11(22) 



T eEx. 



(a^ +/32)cos6'-H2a/3sin6' 

hK^) \ (a^ + l3^) sine + 2aP sin 9 cos 9 



(A4) 



Si[/(k)]l2(21) = ± 



i(a2-/32) 



2hKg 

dfiE+k) df{E^k) 



E ^^""-y 



x,y 



sm U 
- cos( 



dE- 



+k 



dE- 



(A5) 



Using Eqs. (|A4p - (jA5p and the velocity matrix elements 
©-([H]) one can easily prove, that the spin-incoherent so- 
lution (HH) with /12 = /21 = satisfies Eq. (|25|) . 
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